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SUMMARY 



The relations between load on the structure and rota- 
tion of a joint can "be used to estimate the lowest criti- 
cal load after the equation for neutral stability has been 
tested for three assumed critical loads, each of which is 
less than the lowest critical load. 

The solutions of six simple problems are included to 
illustrate the application of the method of estimating 
critical loads and to reveal certain characteristics of 
the method that should be known by the practical engineer 
using it. Four of these problems are concerned with mem- 
bers that lie in the elastic, or long-column, range. The 
other two problems are concerned, with members that lie in 
the short-column range. 



INTRODUCTION 



One of the problems in the design of structures is to 
make certain that the compression members are stable un- 
der the loads to be carried. For structures built with 
the members joined to each other by frictionless pins, the 
usual column formulas can be directly applied to the de- 
sign of the compression members. For structures built 
with the members continuous at the joints, however, the de- 
sign of any one member is dependent upon the design of all 
other members. 

Reference 1 shows how the principles of the Cross 
method of moment distribution can be used "to check the 
stability of structural members under axial load and hence 
the safety of the design. In this method, the critical 
load for the system of members is calculated and compared 
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with the applied load. If the critical load is greater 
than the applied load, the members are stable. If the 
critical load is less than the applied load, the system is 
unstable'" and a larger size for one or more' of the compres- 
sion members must he selected. 

One disadvantage of any method of calculating the crit- 
ical load for a system of structural members under axial 
load is that, for each type of instability, there is a cor- 
responding critical load. In design, the lowest critical 
load is the only one of interest. When the stability of a 
group of structural members is checked, it is therefore the 
lowest of these critical loads that must bo calculated and 
compared with the applied loads. 

Although the two equations for neutral stability given 
in reference 1 are algebraic in appearance, they are fun- 
damentally transcendental in character with the unknown 
critical load entering in angles. The m-ethod of solution 
used in reference 1 was to assume several values of the 
critical ioad and to test one of the equations for neutral 
stability. If the first load in the series of assumed 
loads is made sufficiently small, the lowest assumed load 
that just satisfies this equation is the critical load de- 
sired. 

* 

Unless the designer is fortunate in selecting the as- 
sumed critical loads, considerable time and labor are re- 
quired to find the lowest critical load. In order to make 
the theory of reference 1 more useful in practical calcula- 
tions, a method of estimating the lowest critical load is 
presented in this report. 

In this paper, as in reference 1, it is assumed that 
the members lj.e in a plane and that buckling occurs in this 
plane. It is further assumed that the joints of the struc- 
ture are held rigidly in space but are free to rotate un- 
der the elastic restraint of the interconnecting members* 



DEFINITIONS AND SYMBOLS 



The following definitions of stiffness and carry-over 
factor are the same as those given in references 1 and 2. 
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Stiffness 

If a member is on unyielding supports at each end, 
the moment at one end necessary to produce a rotation of 
one-fourth radian of that end is called the "stiffness." 
The stiffness of a member will depend upon the amount of 
restraint at the far end. In the derivation of the crite- 
rion for stability as given in reference 1, three types of 
restraint at the far end are considered. She symbols used 
to designate the 'stiffness for the different types of re- 
straint are: 

S, far end fixed. 

S 1 , far end elastically restrained. 

S" , far end pinned. 

The stiffness of a member computed according to the 
foregoing definition is one-fourth that computed accord- 
ing to the definition given in references 3 and 4. In the 
Gross method, the relative stiffness of the members is of 
importance and not the absolute value. The foregoing def- 
inition was selected so that the stiffness of a member of 
constant cross section_with no axial load_and fixed at 
the far end would be El/L instead of 4BI/L. 



Garry-Over Factor 

If a member is on unyielding supports at each end and 
a moment is applied at the near end, the ratio of the mo- 
ment developed at the far end to the moment applied at the 
near end is called the "carry-over factor." As in the 
case of stiffness, the carry-over factor will depend upon 
the degree of restraint at the far end of the member; The 
symbols used to designate the carry-over factor for the 
different types of restraint are: 

C, far end fixed. 

C', far end elastically restrained. 



C" = 0, far end pinned. 
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Sign Convention. 



The sign convention used in this report is the same 
as that used in references 1, 2, and 4-. A clockwise mo- 
ment acting on the end of a member is positive. A coun- 
terclockwise moment acting on a joint ie positive. i An " ex- 
ternal moment applied at a joint is considered to act on 
the joint. A positive moment- acting on the. end of a raom^ 
ber causes positive rotation of that endi 



B , modulus of elasticity. 

mm 4 * 

E, effective modulus of elasticity. 

I, moment of inertia of cross section of member 
about a centroidal axis normal to the plane 
of "bending. 

L, length of member. 

W, total load on the structure. 

P, axial load in member (absolute value). 

A, area of cross section* 

c, restraint coefficient in the usual column formula. 



Symbols 
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METHOD OP ESTIMATING- THE CRITICAL LOAD 



The method of estimating the lowest critical load is 
"based upon the principles discussed in references 5 and 6 
.for the analysis of experimental observations in problems 
of elastic stability. In reference 5, Southwell mentions 
that the unavoidable imperfections in practical structures 
prevent the realization of the concept of a critical load 
at which deflections begin. Instead, the initial deflec- 
tions present in practical structures steadily increase 
with load and, according to the usual theory, the deflec- 
tions become infinite as the critical load is approached. 

In references 5 and 6, the relation between load and 
deflection for problems of elastic stability is also dis- 
cussed. The more general -relation given in reference 6 

y - y, 

shows that, if — — — is plotted as ordinate against 

y — y 2 as abscissa, the curve obtained when P approaches 
? crit is essen * ially a straight line the inverse slope of 



which is P 



cr 



it P i ' where 



y is deflection at axial load P in a member. 

y and P , initial values of y and P, respec- 
1 1 tively. 

P crit* lowest critical load. 

P x < P < P crit '" <1> 

Thus, if simultaneous readings of load and deflection re- 
corded in a test are plotted as just described beginning 
with any load F 1 as the initial reading, the value of 

P crit ~ P i is readily obtained. The value of . ? cr it is 
then §iven by the relation 



P 



crit = < P crit - P i> + P i ( 2 ) 



The relation between load and deflection can also be 
applied to load and rotation of a joint. In order to use 
this relation in theoretical calculations, there must be 
initial rotation of the joints. This rotation is obtained 
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by a fictitious external moment M applied at some joint, 
after which the load on the structure ia imagined to "be 
applied. The effect of the tension. in the tension members 
is such as to reduce the rotations caused "by the external 
moment. !£; whereas the effect, of the compression in the 
compression members is such as to increase the rotations. 
As the lowest critical load is approached, the effects of 
compression completely overshadow the effects of tension 
with the result that the rotations become infinite. 

If tho distribution of the total load W on the 
structure does not change as W increases, then the axial 
load in each member is proportional to W. Thus, if 

e - e 

^ ^— is plotted as ordinate against 8 - 8 X as abscissa, 

the curve obtained when approaches ^crit iB osson- 

tially a straight line the inverse slope of -which is 
w crit ~ w i» where 

8 is rotation of a joint under the moment M 
at load TP on the structure. 

& 1 and W x , initial values of 6 and W, respectively. 
^crit» lowest critical load. 



and 



TT X < ff < TT crit (3) 



Thus, if simultaneous values of load and rotation are plot* 
ted as just described beginning with as the initial 

load, the value of ^cr±t " is easily obtained. The 

value of tf cr ifc is then given by the equation 

*crit = (*crit - ff x> + W i < 4 > 

The procedure to be used in estimating tho critical 
load for a group of structural members is: 

1. Assume three values .of W that are known to 
be less than ^ cr if Ekis condition is satisfied if 

the values of W are selected so that the axial load 
in each compression member is less than the strength 
of that member with both ends pinned. 
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2. Calculate the rotation 6 of some joint for 
each of the assumed, loads W "by use of equations 
given later. 

3. Designate the lowest assumed value of 

and the corresponding value of 6 as W 1 and 8^ 
respectively. 

6-6, 

4. Plot the curve of as ordinate 

against 6 - 8 X as abscissa. The three assumed loads 

will give two points on this curve, which are suffi- 
cient to establish an approximate value of ^ C rit ~* 
TT X and, hence, of ^ cr it* * n Poetical calculations, 
the actual plotting of the curve can be omitted be- 
cause the inverse slope ^ cr it ~ "\ would always be 

calculated from the numerical values used to plot the 
curve. If more than three values of W are assumed, 
however, it may be of interest actually to plot the 
curve. 

ROTATION OF A JOINT 



The rotation 8 of a Joint is easily calculated by 
the methods of moment distribution. Either of two equa- 
tions may be used, according to whether the stiffness or 
the series criterion for stability forms the basis. (See 
reference 1 . ) 

Stiffness Criterion for Stability 

Assume that an external moment M is applied at joint 
b in figure 1. The moment -M added to balance this 
joint is distributed among the members as follows: 



to member bc x 



* 5 'be, 
— — j — ~=- 1 to member bc s 
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The moment distribution analysis is now complete as far as 
moments at joint' b are concerned. (See corresponding 
discussion in reference 1. ) 




c l 



\ II, external moment 




Figure 1 



According to the definition of stiffness, the moment 
distributed to any member must "be the rotation of the joint 
multiplied by the stiffness of the member. Hence, 

M_ 



be 



is the rotation in quarter-radians of joint b 



caused by the external moment M. For the purpose of es- 
timating critical loads, M can have. any finite value. 
For the most convenient value, M = - 1, the rotation 6 
is, in quarter-radians, 



6 = 



s s, bc 



(5) 
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Por stability, the moment in the members and the ro- 
tation of the Joints must "be finite. The stiffness cri- 
terion for stability is therefore (reference l) 

S S'^ > 0 (6) 

The condition of neutral stability gives the critical 
"buckling load for the structure and is obtained by setting 
the stiffness stability factor £S »ij C equal to zero, or 

Z S^ c = 0 (?) 

Formulas (6) and (7) are also derived in reference 1. 
Che critical load for the structure is obtained by testing 
equation (7) for different assumed critical loads; the low- 
est assumed critical load that just satisfies equation (7) 
is the critical load desired. If the applied load is less 
than this lowest critical load, the structure is stable; 
if not, the structure is unstable. The method of estimat- 
ing the critical load is therefore a tool to aid in find- 
ing the lowest critical load that will satisfy equation (7). 



Series Criterion for Stability 

Assume that an external moment M is applied at joint 
b in figure 2. Jrom the corresponding analysis in refer- 
ence 1, it follows that the total moment in members' ba at 
joint b is 

U 2 S ' 



ba /, . „ , ,.8 . „3 



or 



S bc + £ S 'ba 
M S S 



(1+ r + r s + r 3 + ....) 



ba 



■where 



S bc + 2 S« ba 1 - 



r _ _£i>c_£bc s ch G o 3a , R s 

S bc + S S»^ a S^. + 2 S« cd <8J 
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Figure 2 



According to the definition of stiffness, the total 
moment in members "ba at joint b must be the rotation 
of joint b multiplied by the total stiffness of members 
ba. Hence, 

M 1 

* s bc + S S 'ba 1 - r 

is the rotation in quarter-radians of joint b caused by 
the external moment M. For the purpose of estimating 
critical loads, M can have any finite value. For the 
most convenient value, M = - 1, the rotation 6 is, in 
quarter-radians , 

1 1 

6 = ■ (9) 

S bc + 2 S 'ba 1 ~ r 

For stability, the moment in the members and the ro- 
tation of the joints must be finite. As stated in refer- 
once 1, the series criterion for stability is therefore 

r < 1 (10) 

The condition of neutral stability gives the critical 
buckling load for the structure and is obtained by setting 
the series stability factor r equal to unity, or 



r = 1 



(11) 
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formulas (10) and (ll) are also derived in reference 
1. These expressions are sometimes more convenient to use 
than the corresponding formulas (6) and (7). In cases 
when the structure is symmetrical about a Joint, the ex- 
pressions concerned with the stiffness criterion usually 
involve fewer calculations. When the structure is symmet- 
rical about a member, the formulas' concerned with the se- 
ries criterion offer certain advantages. Experience in 
the solution of practical problems will dictate which ex- 
pressions result in fewer calculations. In any case, ei- 
ther set is correct and the method of estimating the crit- 
ical load is a tool to aid in finding the lowest critical 
load that will satisfy the equation for neutral stability, 
either equation (7) or equation (11 ). 



In order to check the stability of a group of struc- 
tural members by use of the previously given formulas, 
additional equations for the carry-over factor and the 
stiffness are required. 



GARRY-OVER EACI0R AND STIFENESS 




Figure 3 



12 



N.A.O.A. Technical. Note Np.. 7X7.. 



Consider the member ij ■ shown, in figure 3, simply 
su.pport.ed at i and. elastically restrained at J "by mem- 
bers ■ Jk, The members jk are also elastically restrained 
at their far ends k. By a- moment-distribution analysis 
given in reference 1, the carry-over factor Is 

' (12> 

and the stiffness s 'ij is 

Substitution of equation (12) in (13) gives 

S" 

id , . 

FS^ 



For member ij , the limiting values of the carry-over 
factor and the stiffness given by equations (12) and (14), 
respectively, are obtained as follows. When the far end 
j is pinned, there is no elastic restraint at j and 
2S 'jk = °» ror 'this limiting condition, the carry-over fector 

0 'i j = °"i,J = 0 and the stiffness S = S"^. When the 
far end j is fixed, there is complete restraint at j 
and ES^ic = 00 . Tor this limiting condition, the carry- 
over factor C '.j^ = O^j and the stiffness = S^j , 
where 

S» ^ • 

15 1 - c J-i °ij 



Up to this point, all the equations in this report on 
stability are general* In nearly all of the cases encoun- 
tered in practice, however, the cross section and the axial 
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load do not vary along the length of each member. 3Por 
this special case, = 0j it s " i j = S "ji» an4 S ij = 

Sji. In practical problems, the numerical values for; 

these quantities are obtained by. use of the tables given 
in reference Z. _ 

PROBLEMS 



The purpose of including problems is to demonstrate 
the previously described method of estimating critical 
loads. Sir simple problems have been selected to reveal 
certain characteristics of the method that should be known 
by the practical engineer using it. In order to show the 
accuracy of the estimated critical load, the correct value 
of the critical load for each problem is first established. 

The tables of reference 2 were used in the numerical 
evaluation of the stiffness and the carry-over factor. Al- 
though interpolation in these tables is unnecessary for 
the solution of practical problems, interpolation was used 
for the solution cf these problems to show clearly how the 
estimated critical load becomes more accurate as the as- 
sumed loads W and W x approach W Qrit . 

In problems 1 to 4, it is assumed that the members 
are subjected to low stresses corresponding to the elastic 
range where the effective modulus E is equal to Young's 
modulus E. In problems 5 and 6, the compression members 
are loaded above the elastic range where E < E. In other 
words, for problems 1 to 4, the compression members lie in 
the long-column range ; whereas, in problems 5 ahS""6, the 
compression members lie in the short-column range. 



Problem 1 



Problem: To calculate the critical load for the pin- 
end strut shown in figure' 4. 



Compression P 



EI = 10,000 lb. in. 3 



7\ 



•L = 100 in. 



Figure 4. - Problem 1. 
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The equations concerned with, the stiffness criterion 
for stability are used in the solution of this problem. 
Imagine the external moment M to be applied at Joint b» 
The correct value of the critical load is therefore the 
lowest assume'd load that will satisfy equation (?).. There 
being only one member be, the summation sign is omitted. 
Because this member be is pinned at the far end .s, the 
single prime on S-^ is replaced by a double prime. (Thus 
for this problem, equation (7) becomes 

S» bc =0 (16) 

For member be 

EI = 10,000 lb. in. 3 
L = 100 in. 

Consequently , 




From the tables of reference 2, it is found .that the small- 
est value of P to satisfy equation (16) is the value of 
P giving l/j = tt. Therefore the correct critical load 
is 

P crit = tt S = 9.87 lb. 

which agrees with the value given by the well-known Euler 
column formula 

Pcrit - ^¥ (I?) 

The estimated value of the critical load is given by 
the inverse slope of the approximately straight line ob- 

6-8 

tained by plotting — ■ as ordinate against 8-6, as 

abscissa. Tor this problem, TT = P and equation (5) be- 
comes 



(18) 
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The values of 8 are given in table I for a series of 

6-6 

assumed loads P; and values of 8 - Q x and — are 

given for P x = 0, 3, and 7 pounds. Table I was made ex- 
tensive in order to show how the estimated critical load 
is affected "by P x as well as by the values of P at 
which the inverse slope is computed. 

The approximately straight lines that correspond to 
Pj^ =0, 3, and 7 pounds are plotted in figure 5. Inspec- 
tion shows the lines corresponding to P x = 0 and Pj^ = 
3 pounds to he essentially straight. As only two points 
establish the line for P x = 7 pounds, no conclusion re- 
garding its straightness is justified. 

If P x = 0, then the inverse slope between P = 1 
and P = 2 pounds is (see table I) 

P - P = 0^002202 - 0.000993 m 10 . 33 lb- 

crit i 0.001101 - 0.000983 

from which 

p crit = < P crit - V + P x = 10 * 33 + 0 = 10' 33 . 

The results of a number of calculations of this type 
for other values of P and P x are given in table II. 
Inspection of this table shows that, for any value of Pj , 
the estimated critical lpad becomes more accurate as the 
values of P between which the inverse slope is calcu- 
lated approach Pcrit* ^ke accuracy is also increased as 
P x approaches Pcrit« 

Problems 2, 3. and 4 

The purpose of problems 2, 3, and 4 is to study the 
effect of the tension in tension members on the estimated 
critical load for a structure. In these problems, end b 
of the strut -used in problem 1 is restrained against rota- 
tion, by the adjacent member, ba , which is the same size 
as member be. (See fig. 6.) In problem 2, member ba ... 
has zero axial load. In problems 3 and 4, there is axial 
tension in ba of magnitude P and 3P, respectively; 
In each problem, 'the compression in member be is of mag- 
nitude P. 
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Pro Diem 



Member "ba 



Member "be 



2 
3 
4 



Zoro axial load 
Tension P 
Tension 3P 



Compression P 
Compression P 
Compression. P 



321 = 10,000 lb. in.' 



SI = 10,000 lb. in.' 



A 



~A 
c 



a 



T 



L. = 100 in. 



L = 100 in. 



Figure 6. - Problems 2, 3, and 4. 



Imagine the external momont M to be applied at 
joint b. The correct value of the critical load is the 
lowest assumed load that will satisfy equation (7) which 
becomes, for problems 2, 3, and 4, 

5 "ba + S "bc » 0 (19) 

For each of members ba and be- in problems 2, 3, 
and 4 f 

SI = 10,000.1b. in. a 
L = 100 in. 

On calculation of the values of l/j for each span in each 
of the problems, it ig found by trial that the lowest value 
of P satisfying equation (19), or P cr i t , is 

Probl em .. F crit ( llcr * ) 

2 13.89 

3 15.41 

4 16.93 

The estimated value of-the critical load is given by 
the inverse slope of the approximately straight line ob- 

6 ~ e i 

tained by plotting pT**p"*" as ordinate against 6 - 

as abscissa. for problems 2, 3, and 4, equation (5) be- 
comes ' 
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1 



(20) 



ba 



+ S" 



be 



In tables III, IT, and V, the assumed values of P 
and the corresponding values of 6 are given for problems 
2, 3j and 4, respectively. Che curves established by the 
data in these tables are .plotted in figures 7, 8, and 9, 
respectively. A summary of the corresponding estimated 
critical loads for each problem is given in tahles VI, VII, 
and VIII. 

In figure 7, the curve for P t = 0 is noticeably con- 
cave upward; whereas, in figures 8 and 9, this curve is def- 
initely concave downward. This change from concave upward 
to concave downward is caused by the tension in member ba f 
which results in an o verestimat ion of the critical load 
when the tension in member ba is zero but an underestima- 
tion when the tension is equal to P and 3P. (See tables 
VI, VII, and VIII.) These same conclusions hold in a lesser 
degree when 0 < P l < P orlt . 

When the method of estimating critical loads is applied 
in the solution of practical problems, it is desirable to 
know whether the true critical load is overestimated or un- 
derestimated. ]?rom problems 1 and 2 it is concluded that, 
in the ahsence of tension members, the estimated critical 
loads are all greater than the true critical load. (See 
tables II and VI.) From problems 3 and 4 it is concluded 
that, in the presence of tension members, the estimated 
critical loads are all less than the true critical load. 
(See tables VII and VIII.) The region within which all es- 
timated critical loads are in good agreement with the true 
critical ioad cannot be definitely established in the gen- 
eral case. 

Qualitatively, the region of transition from over- 
estimating to underestimating the critical load can be es- 
tablished by noting the trends in problems 2, 3, and 4. 
In prohlem 2, the poorest estimate of the critical load 
(small values of P and P x ) is 12.. 8 percent on the unsafe 
side. In problems 3 and 4, the poorest -estimates are 69.6 
and 122.7 percent, respectively, on the safe side. In 
problem 2, no tension member is present. In problem 3, the 
size, the axial load, and the number of the tension members 
are the same as for the compression members. In problem 4, 
the axial load in the tension memher is three times the 
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axial load in the compression member. It is therefore 
concluded that the transition from overestimating to under- 
estimating the critical load will occur when the size, the 
axial load, or the number of the tension members is small 
relative to the compression members. 

When all members in a given problem are compression- 
members or when the effects of tension-, are neglected in 
the calculation, as is' sometimes done in stability prob- 
lems of this type, the agreement of the estimated critical 
load with the calculated critical load will be &6 good as 
that found in problems 1 and 2. When the effect of. ten- 
sion in tension members is considered, the precision of 
the ■ estimated critical load can be determined qualitative- 
ly by reference to problems 2, 3, and 4. 

As in the case of problem 1, the agreement of the es- 
timated critical load with the calculated critical load 
for. problems 2, 3, and 4 becomes closer as -the values of 
P between which the inverse slope is calculated, approach 
]?crif The precision also increases as P x approaches 

p crif ( See tables VI, VII, and VIII.) 

TEE EFFECTIVE MODULUS 



Before the theory of this report can be applied to 
problems involving compression members that are stressed 
beyond the elastic range, as in problems 5 and 6, it is 
necessary to introduce an effective modulus E so de- 
signed that the results will be in good agreement with tho 
accepted column formulas. 



Compression MemberB 

Most engineers are familiar with the origin of the 
accepted column curve for a given material. At low 
stresses (stresses less than about one-half the yield 
point of the material )■» the column strength is given by 
the Euler formula. At high stresses, laboratory tests 
always show that the column strength falls short of the 
value given by the. Euler formula. An empirical straight 
line or a parabolic curve is sometimes used to give tho 
column strength within this range. 
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The theory of this report gives a "buckling load that 
is analogous to the strength given "by the Euler column 
formula. As in the case of the Euler formula, these loads 
would not check experimental values at high stresses. A 
reduced strength must therefore "be calculated consistent 
with the accepted column formula for the material of which 
the members are composed. These calculations are best made 
by use of the effective modulus E = TE. 

Consider the case of an ordinary column. If the Euler 
formula is written 

P TT a TE (21) 



A 



<b=k) 



c p 



it will give the strength at both low and high stresses. 
At low stresses, t = 1» whereas, at high stresses, t < 1. 
The problem is to determine how the effective modulus tE 
varies with the stress P/A. - - 

If equation (21) is solved for TE , the following 
equation is obtained 



TT A f- M}/ 



The accepted column formula for any material is always 
given in terms of the effective slenderness ratio 
(L/ V~c" p). Thus, if any one of these formulas is solved 
for ( L/ J c p ) and this value is substituted into equa- 
tion (22), there results an equation for the effective 
modulus tE that is a function of the stress P/A. 

For example, consider the case of S.A.E. 1025 steel. 
The column formulas for this material are: 



p 

Eor j < 18,000 lb. per sq. in. 

■ £ = TT S E 

A "H-) a 

Eor 36,000 >| > 18,000 lb. per sq. in., 



(23) 



? = 36,000 - 1.172 i (|) (24) 



A 
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If equations (23) and (24)' are solved' for (L/ p) and 
these values are substituted' into . equat ion (22) , the fol- 
lowing values' are obtained for the effective modulus E si- 
te .' 

For j- < 18,000 lb. per sq. in., 

I = TE = E (25) 
For 36,000 > | > 18,000 lb. per eq. in., 

/36.000 - f\ 

5 . TE - \ f ( ,^^4: ) (26) 
it A \ 1.172 / 

Equation (25)' shows that, in' the long-column or elas- 
tic range, E = E . Equation (26) shows that, in the short- 
column range, 2 is a function of the stress P/A in the 
member and is in no way dependent upon the stiffness or end 
fixity of the member. 

When the compression members of the structure are 
stressed beyond the elastic range, the methods outlined in 
this report can also be used to calculate the critical 
load. The procedure is the same as in problems 1 to 4 ex- 
cept that, for each assumed load W on the structure^ 
there is a different value of the effective modulus E. 
These values of 1 are obtained by use of equations (25) 
and (26) if the material is S.A.S. 1025 stool. For any 
other material, corresponding equations can bo derived. 



Tension Members 

When the effect of axial load in the tension members 
is considered, the variation of B" with stress for ten- 
sion members can be established, theoretically, by tho use 
of the double-modulus theory of bonding and of tho stress- 
strain curve of the material. Tor such calculations, how- 
ever, the stress-strain curve must be accurately drawn to 
a suitable scale. In the absence of a known or a calcu- 
lated variation of E with stress, the_f ollowing approx- 
imate method can be used to establish E for tension mem- 
bers t 

1. When the stress is v less than the maximum al- 
lowed for a column of the same material, use the same 
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values of 1 for tension as for compression at the 
same stress. 

2. When the stress is greater than the maximum 
allowed for a column of the same material, assume 
' that E = 0 . 

The -values of E for tension members obtained . by 
this method will he conservative. Whether or not they are 
too conservative is a matter to he settled by tests. In 
the regions of yield point and of maximum tensile strength, 
the flatness of the stress-strain curve will certainly 
cause E to approach zero. Because the maximum stress al- 
lowed in columns is closely associated with the yield 
point, this method offers a convenient solution of E" for 
tension members. 



Problems 5 and 6 

The purpose of problems 5 and 6 is to show that the 
method of estimating critical loads presented in this 
paper gives good results when the compression members lie 
within the short-column range. Except for the different 
dimensions and the fact that the members with axial load 
are stressed beyond the elastic range, these problems are 
similar to problems 2 and 3, respectively. 

Problem Member ba Member be 

5 Zero axial load Compression P 

6 Tension P Compression P 



L = 60 in. > 


„ L = 50 in. ^ 


c L = 50 in, > 


< L = 50 in. ^ 







•A 



For -problem 5 



^ For problem 6 



Material: S.A.E. 1025 steel tube continuous from a to 
c with the following dimensions! 

Diameter, d 1.625 in.. 

Wall thickness, t .065 in, , 

Area, A ,3186 sq.. in. 

Moment of inertia, I .0970 in, 4 

Figure 10. - Problems 5 and 6. 
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The essential dimensions for problems 5 and 6 are 
given in figure 10. The effective modulus I for any 
memher is a function of the stress P/A in that member. 
The numerical value of . § for any assumed load P is 
therefore given "by equations (25) and (26), the material 
"being S.A.E. 1025 steel. By the same methods as used in 
the solution of problems 2 and 3, it is found that the 
lowest value of P to satisfy equation (19), or ^crit» 
is 

Problem ?crl A(lb. ) 

5 9,420 

6 9,510 

The necessary calculations for estimating the critical 
loads for problems 5 and 6 wore made by the sane methods 
used for problems 2 and 3 except that in the calculation 
of the stiffness of the members the effective modulus 1 
was used in place of Young T s modulus E. The results of 
those calculations are given in tables IX to XIV, 

6 - 6, 

In figures 11 and 12, - is plotted against 

9 - 6j for problems 5 and 6. It is preferable, however, 
to compare the results given in tables XI and XIV rather 
than to draw conclusions from figures 11 and 12, 

In problem 5, the axial load in member ba is zero; 
whereas, in problem 6, member ba is subjected to axial 
tension equal to the axial compression of member be. 
Comparison of the estimated critical loadB for each of 
these problems (tables XI and XIV) shows that the critical 
load is usually, but not always, overestimated when the 
tension in member ba is zero and is usually, but not al- 
ways, underestimated when the tension in member ba is 
equal to the compression in member be. These same conclu- 
sions wave found for problems 2 and 3. 

Comparison of the precision of the estimated critical 
loads for problems 5 and 6 (tables XI and XIV, respective- 
ly) with the precision of the estimated critical loads for 
problems 2 and 3 (tables VI and VII, respectively) is not 
justified. For problems 2 and 3, the series of estimated 
critical loads are based upon values of P taken at inter- 
vals of roughly 10 percent of ^c^it? whereas, for problems 

5 and 6, this interval was not maintained. When the members 
lie in the short-column range, an estimated critical load 
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based upon small values of P, which lie in the elastic 
range, gives an estimated critical load much higher than 
the true critical load. Consequently, in problems 5 and 
6, the assumed values of P for' which the estimated crit- 
ical load was obtained were made to correspond to values of 
P/A that are associated with the short-uolumn range. 

In the solution of any problem, it is necessary only 
that the assumed loads be less than the true critical load. 
As the assumed loads approach the true critical load, the 
precision of the estimated critical load is increased. 
(See tables XI and XIY.) It is therefore desirable to ex- 
orcise the best judgment possible in the selection' of the 
assumed loads. . In any case, however, the method of esti- 
mating tho critical load as described in this paper '• should 
be regarded as a tool to be used in finding tho lowest 
critical load that will satisfy the equation for neutral 
stability. If it is desired that the estimatod critical 
load be conservative rather than err on the unsafe side, the 
effect of the axial load in the tension members should be 
considered in the calculation. 



1, If the distribution of the total. load W on the 
structure does not change as W increases, then the axial 
load in each member is proportional to W. Thus, if 



the curve obtained when W approaches ^ c rit ■ is essential- 
ly a straight line the inverse slope of which is TC C rit - 
^ , where 

6 is the rotation of a joint under tho moment M 
at load W on the structure. 

Q-l and ^i, initial values of 8 and W, respec- 
tively. 

^crit' lowest critical load. 



CONCLUSIONS 




is plotted as ordinate against 



8-8, as abscissa 



and 
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Thus, if simultaneous values of load and rotation aro plot- 
ted as just described beginning with W3. as the initial 
load, the value of W cr it - ^1 is easily obtained. The 
value of W C rit is then given by the equation 

ffcrit = < ff crit - W i> + ff i 

2. The rotation 0 of a joint can be calculated by 
.the methods of moment distribution. The equation to be 
used depends on -whether the stiffness of the series crite- 
rion for stability forms the basis of the calculation 

3, Por loads within the elastic range, the estimated 
critical load more closely agrees with the calculated crit- 
ical load as the values of W between which the inverse 
slope is calculated approach ^ C rit* ^ e agreement is 
also closor as W^ approaches ^crit* 

4-« 3Tor loads beyond the elastic rango , the results of 
computation have shown that conclusion 3 usually, hut not 
always, applies. In cases where it does not apply, tho 
errors are of the order of a fraction of 1 percent. For 
practical design calculations, . conclusion 3 therefore 
holds for loads beyond the elastic- range as well as for 
loads within the elastic range. 

5, When all members in a given problem are compres- 
sion members or when the effect-s of tension are neglected, 
as is sometimes done in practical calculations, the calcu- 
lated critical load is overestimated. When the effect of 
tension in the tension mombers is considered, however, tho 
calculated critical load is underestimated. The region 
within which all estimated critical loads are in good agree- 
ment with the calculated critical load cannot ho definitely 
established in the general Case. Tho transition from over- 
estimating to underestimating the calculated critical load 
tends to occur, however, when the size, the axial load, or 
the number of tension mombors is small relative to tho com- 
pression members. In many practical problems, tho preci- 
sion with which the estimated critical load. agrees with tho 
calculated critical load can be qualitatively determined 

by reference to tho problems of this report. 

6. The method of estimating the critical load should 
always be regarded as a tool to aid in finding tho lowest 
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load that satisfies the equation for neutral stability. 
This lowest load is the calculated critical load for the 
problem. 

Langley Memorial Aeronautical Laboratory, 

National Advisory Gonnittee for Aeronautics, 
Langloy Field, Va. , June 8, 1939, 
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TABLE. L 



Calculated Dodo, for Problem 1 



p 

lib.) 


© 


P, • Q 

O, a 0. 013333 


f* » » *. 

©,» 0.017010 


P, - 7 lb. 

6, -O.0J3H.S 




e-8, 
P-R 


a-o, 


9-8, 
P-P, 






0 


0.OI3333 


0 


0 










1 


.014114 


.000101 


.000183 










ft 


,015535 


.Qozvn, 


,001)010 










i 


.0I7OTO 


,003757 


.001*5X3 


0 








' 4 


.oiiiss 


.OWJW 


.0014350 


.00*043 


O.OQtOW 






5 




,006103 


.001140* 


.004146 


,ooe+T3 






fa 


.016361 


,01308b 


.OOH747 


,001*71 


.OOSOT3 






7 


,013666 


. 0203S5 


.OOJS1050 


,01457a 


.0O4t44j 


0 




a 


.048467 


.OSS&Ui 


.0044170 


.03J5T1 


.0063158 


.015001 


0./S00| 


i 


.017182 


,08444<j 


.0014054 


.0808% 




,064314- 


.03*157 



TABLE EL 
Calculated Data for FVobltm 2 



(lb.) 



P,- o 

e, « 0.0044447 



9-6. 



p-t; 

v 4-lbJ 



8,» 0.0011434 



0-6. 



p-p, 

\ 41b. I 



P,«9.8 to, 

e, -o.ooiai 



8-8, 



o 

1.4 
6.8 
4v£ 
5.4 
10 
0.4 
*6 



0,0044*67 
.0070077 
.007412? 
•007W4 
.00661)3 
j0O«"S09 
.01014 

.oi3 iai 
.0177/5 
-Oil 91.7 



o 

.MO 34 10 

.00O7J66 

.001*747 

.0011546 

,00*804% 

.0041344 

,0065 14b 

,011046 

■025£b0 



400314157 
AXM.70O7 
.00010318 
.00034104 
,00041203 
.00030648 
,0004447* 
,0001*6*4 
.0060046 



0 

.0004771 
.0016073 
.00*1177 
.903*371 
,0017717 
.063164. 



,00048*11 
00057411 
OOOTOBia 
.100UJW 
,0013160 

.oo&awz 



,0043330 
,010744 



TABLE 32 

Calculated Data for Probkm 3 



.8 



R - o 

ft .0.0066667 



a-e, 



«;», 

P-P, 

jaidlflnf 



• 4.8 lb. 



B-e, 



e-e, 
p-r 

L 4-lb.J 



P, Ml.l ih 
0, » 0.0011641 



0-6, 



e-e, 
P-P, 



0. 0066667 
.0066113 
,0066023 
.0Q61S4I 
.007*102 
.0077664 
,0085470 
.0011641 
,013166 
.QV.5I5 



0 

.0000316 
.0001334 
.000*11*4 
,00062.35 
.0010117 
.0016803 
.0031174 
,0064117 



0.0O0O3O3X 
.0000+075 
.0O0067IH 
,0000 174*3 
,0601*746 
,00011586 
,000*1441 
■00050777 
.0013783 



0 

,0003011 

.0OO177J 

.0015371 

,00*1750 

.0061773 



ojjooissn 

.00014*1) 
,000 SMS 

.00077*14 
.0080331 



0 

.0036023 
.016531 



0,0010014 
.0051720 



table: v: 

Calculated Data for FVobkm 4- 



P, ».o 

6, -0,0066667 



9-9, 
P-P, 

(xadWi 
V. 41b. / 



R »«l lb. 

6, » 0.0056644 



e-e. 



o-o, 

P-P, 



p, » 11.1 fc. 

©, ■ C. Q047O33 



e-e. 



P-P, 



0,0066647 
.0061135 
.0058X4) 
.0US6044 
,0056666 
.00574*0 
,P0hO74l 
.0067093 
.0071834 
.0(3763" 



.0O0J3IZ. 
-.0008436 
-.00018*3 
-.OOOIW 
-.0008847 
•.OOO3S06 
.0000366 
.OOI3I67 
,0070171 



0.00OJMH 
.00064811 
-.00011661 
-.00014676 
-.00QI040B 
..tooDsrta 

.oftOtifMt 



0 

,0000140 
,000017k 
,CC*S1IT 

,ooiO(B7 
.OOSMTO 
,0O0t>flM 



•moomm 

jMOOXATM 
.000074004 
.0O0H-W+ 
.000670+7 
.00071118 



0 

,0016801 
.0070*13 



ftOWWSO 
0060741 



TABLE II 

Summary of Estimated Critical Loads 
for Problem 1 



£ p crlt (theoretical) = 9.87 lb. 3 



p 

1 


Values of P 
"between, which 
slope is calculated 
(lb.) 


I'crit (estimated) 


P crit (estimated) 


(lb.) 


(lb.) 


P crit (theoretical) 




1 and 2 


10.33 


1.047 


0 


4 and 5 


10,10 


1.023 




8 and 9 


9.89 


1.002 


3 


4 and 5 
8 and 9 


10.03 
9.88 


1.016 
1.001 


7 


8 and 9 


9.87 


1.000 



TABLE ~SL 
5urr.mo.ri) of Estimated Critical Loads 
for Probkm Z 



t f^tt (tWificaj) - 13.81 IbJ 



(lb.) 


Valuta of P 
between which 
slope Is calculated" 
(lb,) _. . 


P crtr (pftnate<j 
._ (It). 




0 


1.4 And 
SA and TO 
IU and M 


iS.frl 
t4.7fr 
13.1* 


i.ita 

1.0*3 
1.005 


4. a 


5.1, and T.O 
II.Z and li.fr 


14.54 
13.14 


1.041 
1.004 


i.a 


il. fc and IS.fr 


15,11 


» 1,001 



TABLE 211 
Summary ot Estimate Critical Loads 
for ProbUm 3 



[p flrit (theoretical)* 15.41 lb.] 



p, 

(lb.) 


Values of P 
between which 
slope is calculated 

. (lb-) 


(lb.) 


FJrif (estimated) 
fj^theordfcj 


0 


I.b and 3.2 
fr.4 and a.o 
l£8 and 14.4 


4.U0 
11.61 
JS.33 


1>,304 
.172 
.175 


4.8 


fe.4 and 8.0 
IZJb and 14.4 


13. SO 
15.34 


o.OTfr 
.178 


II.* 


|*.8 and 14.4 


15.41 


1.000 



TABLE 

of Estimated Critical Lands 
for ProbUm 4- 



[P^ (theoretical) - I*.l3 lb.] 



P. 

(lb.) 


Values of P 
b«tw«en which 
slope Is calculated 
(lb.) 


(lb.) 




^.itOburttM 


0 


1.1 and 3.4 
6.8 and 6.5 
13.6 and 15.3 


-3.64 

e.frfr 

15.75 


-o.ttl 
.1ST 
,130 


J.I 


fr.fl and 8.5 
I3.fr and |5J 


fl.OA 
lb. 18 


oA77 
•Wfr 


11.1 


I3.b and 13.3 


lfr.27 


0.1 H 



TABLE 21 



5urnmary of Estimated Critical Loads 
for Problem 5 



[P w , lt 9Werical)-1 1 +zo lb, ] 



P! 

Ob) 


\Axluu of P 
between which 
slope ia calculated 


P^tjcaWted) 
(.lb.) 


SrBjajjggfcdj 




5,7348 and 6372.0 


11,715 


1.273 


0 


1007.2. and 7>4fr.+ 


1,712,0 


1.031 




8/83.(1 and 6.7COS 


1.40U 


I.0I7 


5734,8 


7,001.2 and TjMfr.4 
6£8Afc and 8,110.8 


1,243.2 
%S3i.5 


.IBS 
1.014 


7^4*.. 4 


0,ZB3:fc and ytz6.0 


1,51.5.5 


1.0(5 



TABLE IE. 
5tiffn»9S ValyfcS for Pi-obkm J 



Membw bo 




P 


£. 

A . 


c 




(Up 




o 


a 


158,000,0110 




0 


40,747 


14,711 


S7J4.S 


10,000 


!UJ ) an,<wo 




UtfU 


1:3,54* 


sr.sta 


*pr»xi 




*7, **4,ooa 


W*707 


u*i» 


W,5M 


Jf4,*M 




Tt, 000 


ow> 


4i,m 




14,414 


sa,lW> 


7*44.4 


lAooo 


t+,»fa,ooo 


4B,S3i 


t.Bi*4 


10,413 


44,»S 




2*,0OQ 


(54,478,000 




A.O0OB 


&.0H.5 


Sfc.oofc 


a,Tto.6 


2JB,ooo 


ft, 34J.OO0 




3.4444 


(0,72* 


«3,o4o 



For mambw ba: P-0, ^"D, E. * £8.KliJ*lb, per scVln., 



table: x 



Calculated Data for Frobterrv 5 



p 

ObO 


e 

^radlom^ 


P,«0 

0,-0.13380X10"* 


P, - $734.6 lb. 


FJ < 7^4fc,4 lb. 

0,'Cvsa.MOK.to" 4 


9-9,. 


e-e, 


e-e, 

/ radian' 
\ ? ! 


6-8, 


9-0, 




p-f; 


P-P! 


p-p, 

( radian \ 


(radian) 


TO 


n i 


V 41b. / 


0 

S7S4.0 
M7Z.0 
70W.I 
7444.4 
8MS.4 


.17383 

.10339 

JWJ 

.ztsso 

. 1*7774 

.4^38fl 


0 

.04003xlS 

.04138 

.044*7 

J07I50 

.1431* 


.07731 

-117b* 
JMW 
.33433 


0 

.oirtaontf* 

,0t4H 
.|0»3A 

.24000 


o.«w?»io"' 

.mis 

.24811 . 
.4075* 

■oiwrf 


a 

.was) 


I.U47 



TASLE 331 
Stiffness Values for Problem & 





H«mt>«r3 ba 


and be 




** 


S bft +S bc 


p 


A 


F. 






0 


0 


U, 000,000 


34,«S>1 


o 


40,747 


40,741 


B4.IS6 


3,734.8 


IB, 00° 


23,012,000 


S4,S» 


1.2K.S 


33,404 


13.541 


Ti,744 


^»7a.o 


£0,000 


27,444,000 


54, 707 


1,48 S4 


«,7|*. 


IP.SH 


74, MS 


^OOl.t 


te.ooo 


24,42.7,000 


3-1,413 




33,4SI 


14,44* 


70,047 


J444.4 




24,3? 8,000 


48,834 


3.3124 


JZ.034 


10,413 






16,000 


22,478,000 


43.438 


s.obob 


4f,435 


2,081.3 


31,4*7 


8/120.0 


na.ooo 


IS, 3*6,000 


87.S74. 


3.4444 


45,130 


-10,122. 


34,208 



TABLE SDL 
Calculated Data for Frobkm <» 



p 

Clb.) 


6 


P.-O 

©,»am*44.Rio-* 


P,»J^J4.B lb. 
fl,«0,lt.11fcXl6-* 


P 1 -7 1 fc4fc.4 lb, 

e, »a(toi3 mo"* 


e-e, 


e -0, 


radian j 


ft-a. 


e-e, 


e-e, 


p-p, 


p-f; 


p-f; 

t radial 






e« 






0 

51&4.B 
4372,0 
TJOOT.Z 
74444 

s^w.* 
ajno.a 


.11414 
.»4U> 
.14*7(1 

,iun 

.I733B 
.(1M3 


0 

.OOTSllttf' 
.Oil 1 ) 0 

.utote. 

.05747 
.07004 


OBItfMtlS* 
.0(3817 

,oeswi 

.0470M 
.08SW7 
.170216 


0 

^»t8o 
•05W7 

.U»7 


.10044 
.15705 
.24441 
,SM44 


0 

.ISJUJQ 


(.09735 



TABLE XIV 
Summary of Estimated Critical Leads 
for Prohlem 6 



[ p crit (theoretical) = 9,510 1T>, ] 



(1T>» ) 


Values of P 
"between, which 
slope i.a calculated 
(lb.) 


•^crit (estimated) 
(ID.) 


P crit (estimated) 


■^crit (theoretical) 


0 


5,734.8 and 6 ,372.0 
7,009.2 and 7,646.4 
8,283.6 and 8,920*8 


7,715.2 
8,546.1 
9,442.6 


0.811 
.899 
.993 


5,734*8 


7*009.2 and 7,646.4 
8,283.6 and 8,920.8 


8,761.5 
9 ,532.4 


.921 
1.002 


7,646.4 


8,283.6 and 8,920.8 


9 ,573.6 


1.007 
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Tig. 8 



.007 



.006 
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.0024 



Fig. 9 



.0020 ~ 



.0016 



.0012 

6 -. 9i 
P - Pi 

.0008 



.0004 



*..0004 







































i 


















— 7 

/ 






































r l* 

11.9 


lb. 


















































-A 


/ 


















































D 
















-"-""If 




























/> 

p 




















• 





















.006 



.008 



0 .002 .004 

I'iguxe 9.- Variation of with 6 - 0^ for problem 4. 

(See table V.) 
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gig. 11, 




figure 11.- Variation of ~ " |l with 9 - B 1 for problem 5. 
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